Characterizing correlations with full counting statistics: 
classical Ising and quantum XY spin chains 
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We propose to describe correlations in classical and quantum systems in terms of full counting 
statistics of a suitably chosen discrete observable. The method is illustrated with two exactly solvable 
examples: the classical one-dimensional Ising model and the quantum spin-1/2 XY chain. For the 
one-dimensional Ising model, our method results in a phase diagram with two phases distinguishable 
by the long-distance behavior of the Jordan- Wigner strings. For the quantum XY chain, the method 
reproduces the previously known phase diagram. 



1. Introduction. — Thermodynamic phases are tra- 
ditionally described by correlations of local observables. 
In the simplest case of phase transitions associated with 
symmetry breaking, phases are distinguished by the ex- 
pectation value of a local order parameter. In more subtle 
situations (e.g., Kosterlitz-Thouless phase transition), it 
is the decay of correlations at large distances that dis- 
tinguishes between the phases. In recent years, it was 
realized that other, more sophisticated characteristics of 
correlations may be useful: e.g., the notion of "topo- 
logical order" (involving nonlocal order parameters) or 
entanglement entropy (in the case of quantum systems) . 

In the present work, we consider yet another non- 
local characteristics of correlations based on the full- 
counting-statistics (FCS) approach [J (a related prob- 
lem of order-parameter statistics was studied in Ref. [2]). 
It was pointed out recently (in the context of temporal 
correlations) that analytical properties of the extensive 
part of FCS may be used to distinguish between differ- 
ent thermodynamic phases fT[. Here we apply this idea 
to spatial correlations and illustrate it with two exam- 
ples: the classical Ising and the quantum XY spin chains 
(see also Ref. '?! for an example of one-dimensional free 
fermions, which do not exhibit any phase transition). 

2. FCS characterization of thermodynamic phases. — 
In our approach, a thermodynamic phase is characterized 
by the singularities of the extensive part of a suitably de- 
fined FCS generating function xoW- This construction 
is applicable to any infinite system (either classical or 
quantum, not necessarily one-dimensional) which is peri- 
odic in space and possesses an extensive observable tak- 
ing quantized discrete values (e.g., the number of particles 
or the projection of total spin on a given axis). Consider 
a large subsystem S containing N unit cells of the infi- 
nite system. Let Q be our discrete observable restricted 
to this subsystem and normalized to take integer values. 
Then one can construct the FCS generating function for 
the observable Q, 



where the sum is taken over all integer numbers m and 
Pm is the probability for the observable Q to take the 
value m. The generating function xs(^) ^^^ the form of 
a partition function [5] and therefore must depend expo- 
nentially on the size of the system S [S] : 



Xe(A) (X xo(A) 



N 



(2) 



Here Xo(A) plays the role of the extensive part of xe(A). 
It is periodic in A (with period 27r), but does not have 
to be smooth, and it is the singularities of Xo(A) at real 
values of A that we propose to use as a characteristics of 
the thermodynamic phase [3]. 

Note that the definition above is quite general and ap- 
plies to a vast number of statistical and quantum prob- 
lems (in many situations, one even has a choice between 
different possible observables Q). In Fig. [11 we show two 
such examples: particles on a lattice (with Q being the 
number of particles) and a spin-1/2 chain (with Q being 
the number of up spins) . 

Our FCS construction is related to other characteris- 
tics of correlations. The analytic continuation A — )• ioo 
produces the "emptiness formation probability" (EFF) 
[51 [7]. On the other hand, in quantum systems of non- 
interacting fermions, FCS is known to be related to the 
entanglement [8]. 




Xs(A) = 5]P„,e*^ 



(1) 



FIG. 1; Two examples of FCS in statistical or quantum sys- 
tems. In both examples, the subsystem E is encircled by 
a dashed line, (a) Particles on a lattice, Q is the number 
of particles (in the configuration shown, Q = 5). (b) Spin- 
1/2 chain, Q is the number of up spins (in the configuration 
shown, Q — 2). 
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pretation. If one considers the Jordan- Wigner string 



FIG. 2: (a) The FCS phase diagram of the classical Ising 
chain. The shaded upper region is the "nonanalytic" (NA) 
phase, and the lower region is the "analytic" (A) phase [lOj . 
The phase diagram is symmetric with respect to the change 
of sign of F. (b) The FCS phase diagram of the spin-1/2 XY 
chain at zero temperature. The three phases are: Nonanalytic 
I, Nonanalytic II, and Analytic. On the dashed line 7 = 0, 
\h\ < 1, the dependence of FCS on 7 is nonanalytic. The 
phase diagram is symmetric with respect to the (independent) 
changes of sign of 7 and h. 



Singularities in XoW is a subtle characteristics of 
FCS. In Ref. [3] we argued that they are related to pre- 
exponcntial factors in staggered cumulants of the observ- 
able Q (if the singularity occurs at A = tt). In one- 
dimensional systems, we can interpret hixoi^) as the 
inverse correlation length of the Jordan-Wigner string 
exp{iXQ) (such correlation functions were also discussed 
in the context of integrable systems |S]). In some situa- 
tions, Jordan-Wigner strings may be related to physical 
quantities which are directly observable, e.g., spin corre- 
lations in the example of the XY chain below. 

We propose to classify thermodynamic phases by the 
number and type of singularities of Xo(-'^) at real values 
of A. Thus obtained FCS phases do not have to exactly 
reproduce the conventional phase diagram. Below, we 
illustrate this proposal with two examples. 

3. One- dimensional Ising model. — Our classifica- 
tion results in a nontrivial phase diagram for the one- 
dimensional Ising model. We consider the classical Ising 
chain in an external field described by the Hamiltonian 



H ^jJ2 o-jCTj+i + F ^ cTj , 



(3) 



where the Ising spins aj take values ±1 and the statistical 
weights of spin configurations are given by exp(— i/) (the 
temperature is incorporated in the parameters J and F). 
This model is equivalent to the "weather model" consid- 
ered in Ref. '3', and one finds two different FCS phases: 
"analytic" (A) and "nonanalytic" (NA). In the N A phase, 
Xo(A) has a singularity at A = tt. The phase diagram is 
shown in Fig. [2k, in coordinates F and J [9j [TOj. The 
phase-transition line is given by 



coshF = e 



2, 7 



(4) 



VM) - fl ^fc 



(5) 



(with respect to some reference site Jq), then the FCS 
phase transition corresponds to a nonanalyticity of the 
correlation length of the exponentially decaying correla- 
tion function (T4(0)F7r(j)), as a function of the param- 
eters J and r. In the A phase, the correlation function 
(T4(0)K-(j)) exhibits a pure exponential (at F < 0) or 
a staggered-exponential (at F > 0) decay as a function 
of j, while in the NA phase there are additional incom- 
mensurate oscillations in j. Note that this "FCS phase 
transition" is not a phase transition in the usual sense, 
since the partition function of the Ising model does not 
have a singularity at the transition line. 

4. Spin-1/2 XY chain. — We now turn to a quantum 
example where a nontrivial FCS phase diagram may be 
explicitly constructed: the spin-1/2 XY chain in a trans- 
verse magnetic field. The Hamiltonian of the system is 

m 



H = 



E 



1+7 
2 



1-7 



-J-J+i-K 



(6) 



Without loss of generality, we assume 7 > and h > 0. 
We are interested in the FCS phase diagram with respect 
to the number of up spins (as the observable Q in our 
construction) at zero temperature (in the particular case 
of 7 = 1, the FCS in this system was studied in Ref. T). 
By the Jordan-Wigner transformation, this model can 
be mapped onto a quadratic fermionic system |11L 112) , 
and the generating function xs(A) for a subchain of N 
sites can be written as a A^ x A^ Toeplitz determinant. 
By a simple extension of the derivation in Ref. 7, we find 



Xs(A) = det 

l<j<A:<Ar 



dq 

2^ 



a(g,A)e'«(^-'=) (7) 



with the symbol a{q, A) of the Toeplitz determinant given 
by [cf. Eq. (17) of the first paper of Ref. 1^ 



CT(g,A) = 



1 



1-e*^ 
2 



cos q — h + i^ sin q 

(cos q — /i)2 + 72 sin'^ q 

(8) 

The exponential asymptotic dependence of such de- 
terminants on N is given by the Szego formula |13) . It 
immediately produces the result for xo(-^): 



Xo(A) =exp / -- \na{q,X) . 
Jo ^TT 



(9) 



This FCS phase transition has a simple physical inter- 



provided that a{q, A) has zero winding of the complex 
phase as q varies from to 2tt. In the integral ([9]), the 
branch of the logarithm is chosen by the analytic con- 
tinuation along the real axis of q, and the zero-winding 



condition implies that lncr(27r, A) = lncr(0, A) under such 
an analytic continuation. 

For some values of 7, h, and A, however, the symbol 
([8]) may have winding number one [so that In cr(27r, A) = 
In (7(0, A) + 27ri]. In this case, a modification of the Szego 
formula applies [Tl] : 



Xo(A) 



exp 



27r 



dq 

2^ 



In [cr{q, A) ( 



iqa 



(10) 



where qq is the location of the singularity of a{q, A) in the 
upper half plane of q with the smallest imaginary part. 

A tedious, but straightforward application of Eqs. ^ 
and (10 1 allows us to calculate explicitly xq{X) at all val- 
ues of 7 and h. As a result, we find three FCS phases 
shown in Fig.jSja: two nonanalytic phases and an analytic 
one. Note that the same phase diagram appeared previ- 
ously in the analysis of spin correlations [12 and of EFF 
[7]. The generating functions Xo('^) at typical points in 
each of the three phases are shown in Fig. |3] Below we 
summarize some properties of these phases in terms of 
FCS. 

Nonanalytic I phase (NA I): 7^ + /i^ < 1. In this 
phase, we find (assuming A G [— tt, vr], 7 > 0, ft- > 0) for 
the absolute value of XoW- 



lxo(A)| = 



1 + 7 cos A 
1+7 



(11) 



Expressions for the phase of xo(A) following from Eqs. ([9]) 
and ( 10 ) are lengthy for all the three phases, and we do 



not present them here, except in several particular cases 
where they can be considerably simplified. In the NA I 
phase, the only singularity of XoW is a phase jump at 
A = 7r: 



Im lnxo(i''r) — ± \ tt ~ arccos 



yr 



T 



(12) 



In several special cases, XoW takes a particularly sim- 
ple form. At 7 = 0, 



Xo(A) =exp 



1 



iA I 1 arccos h 

IT 



(13) 



[in this case, the spin chain is equivalent to free fermions 
with the density 1 — (I/tt) arccos h]. 
At h^O, 



Xo(A) 



J^/ 



2 / 1 + 7 cos A 



V 1+7 
At the phase boundary 7^ + ft^ = 1, 

1 



XoiX)=pe'^ + l-p, 



P = ^n 



1-7 

1 + 7 



(14) 



, (15) 



1 Xo(^) 1 




arg x^{l.) - X 


^"V. ^ 













■-■ A \ 


-O.lJi 


\v 


NA II \ 




\^ 


_ NAI \^ 




N 



0.5 



FIG. 3: The absolute value and the phase of the generating 
function Xo(A) at typical points in each of the three FCS 
phases of the XY chain. For better visualization, the linear 
part A is subtracted from the phase of Xo(A). In each panel, 
the three curves are (from bottom to top): NA I phase (solid 
line), 7 = 0.6, h = 0.7; NA II phase (dashed line), 7 = 0.5, 
h — 0.95; A phase (dash-dotted line), 7 — 0.5, h = 1.01. 



which corresponds to independent spins with the proba- 
bility p of pointing up [T^ . 



Note that Eq. (Ill is not even in 7, and therefore Xo(A) 
depends nonanalytically on 7 across the line 7 = 0. 

Nonanalytic II phase (NA II): 7^ + ft^ > 1, \h\ < 1. In 
this phase (again assuming A € [— 7r,7r], 7 > 0, h > 0), 
Xo(A) has phase jumps at points ±Ac given by 



cosXr 



1 



Zl 



h+^j^ + h^-1 



7(1 + Z2)' - 1+^ 

For the absolute value of Xo(A), we find 
lxo(A)| = 




(16) 



(17) 



Once again, we do not present here full expressions for 
the phase of Xo(A). The phase jump at Ac is given by 



Im In 



Xo(Ac + 0) 
Xo(Ac-O) 



= arccos 



2zi 



(18) 



(this jump tends to zero at the phase boundary). At the 
boundary with the NA I phase, Ac — ^ tt, and Xo(A) is 
given by Eq. (15). At the boundary with the Analytic 
phase h = 1, one finds Ac — > 7r/2. In the whole NA II 
phase, Im lnxo(i7r) = ±7r, which implies that Xo(A) is 
smooth at A = tt. 

Analytic phase (A): \h\ > 1. In this phase, Xo(A) has 
no singularities in A. For its absolute value we find: 



lxo(A)| = 



lh+ ^72 cos2 A + ft,2 _ 1 

h + ^J^ + h^^ 



(19) 



Throughout this phase, Im lnxo(i''') = iTr (just like in 
the NA II phase). At 7 = 0, the generating function takes 



the particularly simple form Xo('^) = i^ [a- completely 
filled band of free fermions] . 

We notice a remarkable property of the A phase: 
throughout this phase, XoCtt) = — 1. This means that the 
correlations of the Jordan- Wigner operators ([5]) [where 
tjfc denote now the z components a^] decay slower than 
exponentially. This property was described in Rcf. 1 
as "confinement of dual domain walls" , and it can also 
be understood in the fermionic language. At \h\ > 1, 
fermions form a completely filled (or completely empty) 
band, and the anisotropy terms (involving 7) introduce 
local Cooper pairs, which, however, change the num- 
ber of particles by two. Thus the expectation value of 
the Jordan- Wigner parity string (T4(0)F7r(j)) is only af- 
fected by Cooper pairs intersecting one of the end points 
of the string (0 or j). But, since pairs are local, the num- 
ber of such pairs remains finite for long strings, which 
leads to a saturation of the correlations {V^{0)VTr{j)) at 
large j. 

The PCS phase diagram of the XY chain (Fig. (2)3) 
coincides with that obtained from spin correlations ;12j . 
since, by the Jordan-Wigner transformation, transverse 
spin operators a, are represented by the product of the 
string operator ^ and a fermion operator [11] . As a con- 
sequence, transverse spin correlations are given by the 
Toeplitz determinants which differ from Eqs. ([7| and ([8]) 
only by a shift of the winding number (and by fixing 
A = tt) [12]. Therefore, spin correlations produce the 
same phase diagram, where phases differ from each other 
by the presence or absence of incommensurate oscilla- 
tions in (a^'a^) and by pre-exponential factors in their j 
dependence [H] (the latter is disregarded in our present 
work [S]). 

The phase diagram in Fig. [2]d also resembles those 
based on EFP in Ref. [7| and on entanglement entropy in 
Ref. dHl Indeed, the EFP is given by the same Toeplitz 
determinant ([t]), ^ with A -^ ioo [7], and the entan- 
glement entropy depends on the spectrum of a closely 
related block Toeplitz matrix \W . Therefore the phase 
boundaries which are determined by the geometry of the 
square-root branching points in Eq. (|8|) coincide in all the 
three problems. However the PCS classification contains 
additional details related to the positions of A-dependent 
logarithmic branching points in Eq. (pj). 

Finally, we remark that, even though the XY spin 
chain considered in our work maps onto a quadratic 
fermionic system, it does not obey the theorem on factor- 
ization of PCS for noninteracting fermions of Ref. \Wi the 
phase NA II with a singularity at an intermediate value 
of A would not be allowed by that theorem. The reason 
for this discrepancy is that the corresponding fermionic 
system contains pairing terms |llj . and the theorem of 
Ref. (Hjis not, in general, valid for quadratic Hamiltoni- 
ans with pairing (see also discussion in the supplementary 
material of the last paper of Ref. |8j) . 



5. Conclusion. — We have proposed a classification 
scheme of thermodynamic phases in terms of the analyt- 
ical properties of the extensive part of PCS for a suitably 
chosen discrete observable. This proposal was illustrated 
with two examples where PCS could be calculated analyt- 
ically. In more complicated systems, it may be possible 
to access PCS numerically 

In the present paper, we focused on analytic properties 
of the generating function Xo(A) at real values of A. It 
may also be instructive to analyze, more generally, the 
structure of singularities of Xoi^) in the complex plane 
of e*^. This would provide a connection to the theory of 
EFP [5] [7] and to the analysis of PCS for noninteracting 
fermions in Ref. 1161 
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